Non-Gaussian dephasing in flux qubits due to l//-noise 
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Recent experime nts by F. Yoshihar a et al. [Phys. Rev. Lett. 97, 167001 (2006)] and by K. 
Kakuyanagi et al. ( jcond-m at /0609564 ) provided information on decoherence of the echo signal in 
Josephson-junction flux qubits at various bias conditions. These results were interpreted assuming a 
Gaussian model for the decoherence due to 1// noise. Here we revisit this problem on the basis of the 
exactly solvable spin-fluctuator model reproducing detailed properties of the 1/f noise interacting 
with a qubit. We consider the time dependence of the echo signal and conclude that the results 
based on the Gaussian assumption need essential reconsideration. 
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Two remarkable worksi^ have appeared recently, re- 
porting on measurements of the two-pulse echo signal in 
Joscphson junction flux qubits revealing its dependence 
of the time interval between the pulses. The authors in- 
terpreted their results in terms of the recent theory^ of 
the decoherence in qubits caused by 1// noise and ob- 
tained the qubit dephasing rate, T v , based on the postu- 
lation of the Gaussian statistics of the noise. Although 
this looks like a natural starting point, the importance 
of the echo signal data for understanding the underlying 
mechanisms of the qubits decoherence, calls for careful 
examination of the assumptions built into the theoreti- 
cal description. 

In this paper we develop a theory of the time depen- 
dence of the echo signal making use of an exactly solvable 
but experimentally realistic model. We demonstrate that 
in many practical realizations of l/f noise the results 
based on the Gaussian assumption need to be signifi- 
cantly corrected. We show that deviation of noise statis- 
tics from the Gaussian changes significantly the time de- 
pendence of the echo signal. Using the exactly solvable 
non-Gaussian spin-fluctuator model for the l/f noise we 
analyze the dependence of the echo signal on both time 
and the qubit working point, and compare the obtained 
results with those derived within the Gaussian approxi- 
mation. 

Let us start with a brief review the procedure used 
in Refs. [l]0. In order to separate the relaxation due to 
direct transitions between the energy levels of the qubit, 
characterized by the time T\ , the echo signal is expressed 
as e 



-t/2Ti 



p(t), where t is the arrival time of the echo 
signal. The factor p{t) describes the pure dephasing. 

If dephasing is induced by the noise of the magnetic 
flux in the SQUID loop, the echo signal is determined 
by the fluctuating part of the magnetic flux through the 
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The strength of the coupling between the qubit and the 
noise source, depends on the working point of the 
qubit ^ 

Assuming that the statistics of the fluctuations of 
<I>(i) is Gaussian one can express the decoherence 
rate through the magnetic noise spectrum, S$(w) = 
(1/tt) J °° dt cos cot ($(f)$(0)): 



p(t) = exp 



vli 2 



duj Sg,(u>) 



sin 4 (wt/4) 
(tot/4) 2 



(2) 



This expression is common in the theory of spin reso- 
nance and allows one to find the decoherence rate from 
the noise spectrum once the coupling, is known. For 
the l/f noise spectrum 



S$(u>) = A$/uj , 



one obtains*!* 4 ^ 

p(t) = , 



-(rjt) 2 



| \jA<s, In 2 . 



The experimental results of Ref. were fitted to 

P (t) = ex p[-t/2T 1 -r°A-(r*t) 2 ] , 



(3) 
(4) 
(5) 



and T\ and I ± were extracted from the analysis of the 
echo near the optimal point. The dephasing rate, F*, 
was in its turn found from the fitting and displayed as a 
function of the working point. 

The deficiencies of this procedure become clear when 
we compare the approximate expression Eq. ([4]) with ex- 
act equations that follow from the spin-fluctuator model 
which we will now introduce. 
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Spin-fluctuator (SF) model for 1/ / '-noise. One of 
the most common sources of the low-frequency noise is 
the rearrangement of dynamic two-state defects, fluctu- 
ators, see, e. g., book [6| and references therein. Ran- 
dom switching of a fluctuator between its two metastable 
states (1 and 2) produces a random telegraph noise. The 
process is characterized by the switching rates 712 and 
721 for the transitions 1 — > 2 and 2 — > 1. Only the fluctu- 
ators with the energy splitting E less than temperature, 
T, contribute to the dephasing since the fluctuators with 
large level splittings are frozen in their ground states. As 
long as E < T the rates 712 and 721 are of the same 
order of magnitude, and without loss of generality one 
can assume that 712 ~ 721 = 7/2. The random tele- 
graph process, is defined as switching between the 
values ±1/2 at random times, the probability to make 
n transitions during the time r being given by the Pois- 
son distribution. The correlation function of the random 
telegraph processes is (£(t)£(f + t)) = e~ 7 ' T '/4 and the 
corresponding contribution to the noise spectrum is a 
Lorentzian, j/4ir(u) 2 + 7 2 ). Accordingly, should there 
be many fluctuators coupled to the given qubit via con- 
stants Vi and having the switching rates 7^, the dephas- 
ing of the noise power spectrum is expressed through the 
sum J2i v ili/( uj2 + if)- If tne number of effective fluc- 
tuators is sufficiently large, the sum over the fluctuators 
transforms into an integral over v and 7 weighted by 
the distribution function "P(u,7). Upon the assumption 
that the coupling constants and switching rates are un- 
corrected, the distribution density factorizes, V(v,j) = 
'Pv(u)'P'y('l), and the noise spectral function due to fluc- 
tuators reduces to (u 2 ) S(u>), where 
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uj 2 + 7 2 



(6) 



The distribution, ^(7), is determined by the details 
of the interaction between the fluctuators and their en- 
vironment, which causes their switchings. A fluctuator 
viewed as a two-level tunneling system, is characterized 
by two parameters - diagonal splitting, A, and tunneling 
coupling, A, the distance between the energy levels being 
E = \/A 2 + A 2 (we follow notations of Ref. 0] where the 
model is described in detail) . The environment is usually 
modeled as a boson bath, which can represent not only 
the phonon field, but, e. g., electron-hole pairs in the 
conducting part of the system. The external degrees of 
freedom are coupled with the fluctuator via modulation 
of A and A, modulation of the diagonal splitting A being 
most important. Under this assumption the fluctuator- 
environment interaction Hamiltonian acquires the form 



'A A 

/IF-env — C \ —CT Z - — (T a 



(7) 



where c is an operator depending on the concrete interac- 
tion mechanism. Accordingly, the factor (A/E) 2 appears 
in the inter-level transition rate: 



Here the quantity 70 has a meaning of the maximal re- 
laxation rate for fluctuators with the given energy split- 
ting, E. The coupling, A, depends exponentially on 
the smoothly distributed tunneling action, leading to the 
A _1 -like distribution of A, and consequently V 7 cx 7 _1 . 

Since only the fluctuators with E < T arc important 
and temperatures we are interested in are low as com- 
pared to the relevant energy scale, it is natural to assume 
the distribution of E to be almost constant. Denoting the 
corresponding density of states in the energy space as Pq 
we arrive at the distribution of the relaxation rates as^ 



Vy(l) 



P Q T 



0(7o - 7) 



(9) 



The product PqT determines the amplitude of the 1/f 
noise. Indeed, the integral J V-yd) dsy is nothing but the 
total number of thermally excited fluctuators, At- Con- 
sequently, 



P T = N T /C , C = ln(7o/7mi„) 



(10) 



where 7 m ; n is the minimal relaxation rate of the fluctu- 
ators. In the following we will assume that the number 
of fluctuators is large, so that PqT > 1 or A/t 3> 2ir£. 
After that the Eq. © yields: 



S(u) = - x 

Ul 



1 



27o j-KUl 



OJ < 70 
w > 70 ' 



A^^PoT. (11) 



We see that the SF model reproduces the 1/f noise power 
spectrum ^ for lo <C 7o- The crossover from ui~ l to 
ui~ 2 behavior at u> ~ 70 follows from the existence of a 
maximal switching rate 70. Below we will see that this 
crossover modifies the time dependence of the echo signal 
at times t "/q 1 . 

The SF model has previously been used for descrip- 
tion of effects of noise in various systems ^ 10 ' 1 ^ 12 ' 13 ^ 4 
and was recently applied to analysis of decoherence in 
charge qubit s^jm^j^,^^ Quantum aspects of 
the model were addressed in Ref. |23|. These studies, 
demonstrated, in particular, that the SF model is suit- 
able for the study of non-Gaussian effects and that these 
may be essential in certain situations. 

Now we are ready to analyze consequences of the up- 
per cutoff and the effect of non-Gaussian noise, and 
through this identify the validity region for the prediction 
of Eq. ©. 

In the following we will express the fluctuation of the 
magnetic flux as a sum of the contributions of the statis- 
tically independent fluctuators, $(t) = X^ki£i(t), where 
bi are partial amplitudes while £j (t) are random telegraph 
processes. Consequently, we express the product 
as v 2 A where v = \J (v^b 2 ) . In other words, we include 
the amplitude of the magnetic noise in the effective cou- 
pling constant. 

a. Echo signal in the Gaussian approximation. Sub- 
stitution of the distribution © into Eq. © yields for 
AC(t) = -hip(i) 



1 {E 1 K) = {K/E) 2 l0 {E). 



(8) 



fC g (t) 
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lot < 1 
lot > 1 . 



(12) 
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The subscript g means that this result is obtained in the 
Gaussian approximation. One sees immediately that the 
crossover between the w _1 and uj~ 2 behaviors in the noise 
spectrum does not affect the echo signal at long times 
t 3> Jq . However, at small times the decay decrement 
acquires an extra factor jot, which is nothing but the 
probability for a typical fluctuator to change its state 
during the time t. As a result, even in the the Gaus- 
sian approximation at small times K g (t) oc £ 3 (replacing 
JC g {t) oc t 2 behavior). 

b. Non- Gaussian theory: The echo signal given by 
Eq. |T]) can be calculated exactly using the method of 
stochastic differential equations, if the fluctuating quan- 
tity $ is a single random telegraph process, see Ref. 
and references therein. The method was developed in 
the contexts of spin resonance^ spectral diffusion in 
glasses , 25 ' 26 i 27 ' 28 and single molecular spectroscopy in 
disordered media . 29 i 30 i 31 i 32 Averaging in Eq. (p} is per- 
formed over random realizations of $ and its initial states 
and reflects the conventional experimental procedure^ 2 - 
where the observable signal is accumulated over numer- 
ous repetitions of the same sequence of inputs. 

For a single fluctuator with switching rate 7/2 and 
coupling v we nnd 2 ^ 



pi(*) 



-lt/2 



(1 + p)e^ 2 + (1 - /i)e-^*/ 2 
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where p = \Jl — (v/j) 2 . In the appropriate limits this 
can be expanded to give 



lnpi(t\v,j) 



w 2 7t 3 /48 
lt/2 
v 2 t/2j 



t- 1 > 7,w 
7,£ -1 < v 
7 > v,^ 1 . 



(13) 



Note that the last limiting case here is similar to the 
motional narrowing of spectral lines well know in physics 
of spin resonance.— 

In the case of many fluctuators producing 1 //-like 
noise, A/r 3> 1, the sum of the contributions from in- 
dividual fluctuators gives the echo decay decrement 

K si (t) = - fdvV v (v) { ^7^7(7) Inpi(t|«,7) ■ (14) 



Let us assume that the distribution of v is a sharp func- 
tion centered at some value v. This reduces integration 
over v to mere replacing v — > v in the expressions for 
p(v,t\j). This approximation is valid as long as (v 2 ^ is 
finite. This is seemingly the case, e. g., for magnetic 
noise induced by tunneling of vortices between different 
pinning centers within the SQUID loop. The case of di- 
vergent (w 2 ) is considered in Refs. Using then 
@ for the distribution function and using the appro- 
priate terms from Eq. (TIB")) in Eq. (jT5J) we find the time 
dependence of the logarithm of the echo signal, K s {(t). 
For v <C 70 
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(15) 
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FIG. 1: Dephasing component of the echo measurements re- 
plotted from fig 4a of Ref. Q] away from the optimal point. 
The curves show fits to the SF model, eqs. (|15[1 and (|16|l and 
to the p — e~ r ' law. The fitting took into account all data 
points including those that fall outside the range of the plot 
(e. g. those with p > 1). 



where a w 6. 

At small times t <JC v we arrive at the same result 
as in the properly treated Gaussian approach, Eq (|12jl . 
However, at large times, t 3> v~ l . the exact result dra- 
matically differs from the prediction of the Gaussian ap- 
proximation. To understand the reason, notice that for 
■P 7 (7) oc I/7, the decoherence is dominated by fluctua- 
tors with 7 w v. The physical reason for that is clear: 
very "slow" fluctuators produce slow varying fields, which 
are effectively refocused in course of the echo experiment, 
while the influence of too "fast" fluctuators is reduced 
due to the effect of motional narrowing. As shown in 
Ref. HO, only the fluctuators with v <C 7 produce Gaus- 
sian noise. Consequently, the noise in this case is essen- 
tially non-Gaussian. Only at short times f « f 1 when 
these most important fluctuators did not yet have time 
to switch, and only the faster fluctuators contribute, is 
the Gaussian approximation valid. 

For v 2> 70 we find 



IC si (t) 



j Av 2 t 3 /6 
4j At 



t < v 
t > v 



-1 
1 



(16) 



In this case all fluctuators have v ^> 7, hence the result at 
long times again differs significantly from the Gaussian 
result Eq. (fT2|) . 

Discussion: Here we apply the results obtained above 
to analyze quantitatively the decoherence of the flux 
qubit. Fig. Q] shows the time dependence of the echo 
signal measured in Ref. [l| and a fit based on the SF 
model. We consider first the case of v ^> 70 where the 
SF model predicts a crossover from t 3 to t dependence, 
Eq. (flfj)) . We replace the exact result by the interpo- 
lation formula, K, s{ (t) = j At 3 /(6v- 2 +t 2 /4). Figure Q] 
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also shows the commonly used fit p = e~ r ' , which is 
represented by a straight line with slope 2. The fit to 
the SF model seems slightly better, or at least equally 
good. From this fit we can extract the average change in 
the qubit energy splitting Eqi due to a flip of one fluc- 
tuator, v « 4//s~ 1 . It allows us to evaluate the change 
of flux in the qubit loop induced by a fluctuator flip, 
b = v/v$, since v$ = (l/h)dEoi/d& was measured in 
Ref. [y. Using the experimental values for all parameters 
we get b w 4.2 x 10 _6 $ , where $ is the flux quantum 
while the deviation from the optimal working point was 
« 10~ 3 $o- Hence the flip of one fluctuator changes the 
qubit energy splitting by 0.4%. The fit also determines 
the value of the product joA and within our assump- 
tion 70 < v we find a lower estimate for the flux noise 
amplitude y/A^/^a > 1.3 x 10" 6 . 

Similarly, if v <C 70 we can fit to Eq. (|15p . using 
the interpolation formula )C s {(t) = v 2 At 3 / (Q^q 1 +t/ ln.2). 
We do not include here the JC s f oc t behavior at large t 
since it corresponds to K a { > A> 1. The fitting gives 
70 ~ 30 ^s" 1 , which is then also an upper estimate for 
v. We also get v^iV$o = 1-05 x 10~ 6 . 

The quality of the data does not allow us to deter- 
mine which of the cases v 3> 70 or v <C 70 is realized 
in the experiments. However, in both cases the value 
of A® is similar and close to the value obtained in the 
Gaussian approach pi Besides, even without determining 
which case is realized we obtain that both v and 70 can 



not be larger than 30 fis . This gives an upper estimate 
for the change of flux in the qubit loop as one fluctuator 
flips, b < 3 x 10~ 5 $o- Since 70 is expected to grow with 
temperature, we believe that it would be instructive to 
performe similar analysis of the echo decay at different 
temperatures. 

Conclusions: By introducing the spin fluctuator 
model for 1// noise in the qubit level splitting we have 
determined the time dependence of the echo signal. We 
show that the standard quadratic time dependence in the 
Gaussian approximation Eq. has a limited range of 
applicability, and t 3 or t dependencies are found beyond 
this range. Fitting to the SF model also allows us to 
determine the strength of individual fluctuators, and for 
the flux qubits reported in Ref. [l| the change of flux in 
the qubit loop due to the flip of one fluctuator was found 
to be b < 3 x 10- 5 $ - 
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